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Constitutive Modeling and Thermoviscoplasticity 

Development and solution of coupled thermomechanical equations at elevated 
temperature and/or high strain rates are discussed. Three main considerations 
are presented: development of the coupled thermomechanical equations by means 
of the rational theory of thermodynamics, development of a thermoviscoplastic 
constitutive equation which is congruous with the developed coupled equations, 
and the applicability of the developed equations to the treatment by the 
finite element method. 
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Introduction 


It is well established that finite deformations of a solid body particularly 
at elevated temperature and/or high strain rates represent coupled 
thermomechanical processes, which require the simultaneous solution of the 
coupled balance of momenta and energy equations. A proper development and 
solution of such thermomechanical problems requires: 1) adoption of the 
rational theory of thermodynamics, 2) a comprehensive viscoplastic 
constitutive equation which accounts for the strain rate, temperature and 
hardening effects, and 3) compatibility with the available numerical tools, 
particularly the finite element method. These requirements taken together have 
not been used extensively by researchers in dealing with the coupled 
thermomechanical problems. However, because of the need for stringent accuracy 
when solving practical thermomechanical problems such as in rockets and in 
nuclear reactors, the importance of these requirements is being recognized. 
Inoue and Nagaki [1] and Allen [2] developed coupled thermomechanical 
equations with limited applications to one dimensional problems. Ghoneim [3] 
presented a coupled equations, without hardening effects, and applied them to 
a two dimensional axisymmetric problem of compression of a constrained-ends 
cylinder. Lehmann [4] presented a comprehensive analysis of the development of 
the coupled equations with application to the necking problem in a specimen 
subjected to the tensile test. However, a more realistic constitutive law 
which includes strain rate and temperature effects is needed. 

In this paper development and solution of the coupled thermomechanical 
problems is considered, based on the three requirements listed earlier: The 
coupled thermomechanical equations are developed on the basis of the rational 
theory of thermodynamics; a viscoplastic constitute equation which accounts 
for the temperature, strain rate, and hardening effect is proposed; and the 
computability of the developed coupled thermomechanical equation with the 
finite element method is discussed. 
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Development of the Coupled Thermomechanical Equations 


We assume that the state of a material point is oompletely determined by the 

knowledge of the elastic strain tensor £ e , the inelastic strain tensor 
the absolute temperature T, the temperature gradient VT, and a set of 
internal state variables a 1 , i,® l,...,p. Consequently, the following 
constitutive relations may be postulated: 
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where § is the second Piola-Kirchhoff stress tensor, ^ stands for the 
Helmholtz free energy, s means the specific entropy, and q is the heat flux 
per unit area. 

Upon invoking the axiom of admissibility (i.e., the compatibility of the 
assumed constitutive relations with the fundamental equations of mechanics) 
and when adopting the separabilty of the total Green-Lagrange strain energy 

E = E e + E 1 + YT 
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where Y is the coefficient of the thermal expansion tensor, it follows: 

1) * - ip (£ e , E 1 , T, a 1 ) 

s = : l - — * s (E e , E 1 , T, a 1 ) 
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At this stage, we may postulate 


...( 2 ) 

...(3) 


= P^o + J E e : D 1 * : E e + p c y T (1 - In T) + £ : a ...(4) 

4 

where D is the fourth order elasticity tensor, p is the density at the 
reference configuration, and c v stands for the specific heat at constant 
deformation. The tensor £ is the material property tensor. Substituting (4) 
and the Fourier's law, q = - kVT into equations (2) and (3), we obtain 

-k7 ! T + pc v t = -p( 1 : 0':£ e )T.5:! l .p|:j ...(5) 

and 

jL : j/ + P = §:ct + k > 0 

T 

Equation (5) is the coupled heat equation which together with the balance of 
linear momentum equation constitutes the set of coupled thermomechanical 
equations. It might be worth pointing out that the right-hand side of the 
equation (5) represents the mechanical energy generation; the first term 
stands for the reversible part, the second for the dissipated irreversible 
part, and the third for the stored irreversible part due to microstructural 
effects. 
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The Viscoplactic Constitutive Equation 

Solution of the coupled thermomechanical equations requires closed-form 
expressions for (1.5) and (1.6). These may be taken from any of the state 
variable theories [5-8]. Inhere, a constitutive equation which may be regarded 
as a modification of th Bodner-Partom's power law is proposed. The proposed 
equations are thought to be simple, lucid, and consequently very practical. As 
in the case of the internal state variables theory, the proposed constitutive 
equations utilize two state variables: a kinematic hardening state variable 
which accounts for the "rest" stress, and an isotropic hardening state 
variable which accounts for the "drag" stress. Only the isotropic hardening 
state variable will be considered in this paper. 


After adopting the flow rule, we can show that 
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where is the effective stress ( t E * J j SJ j j )» ^ is the 

deviatoric stress tensor, n is a strain rate sensitivity parameter, and Y is 
the isotropic hardening state variable which is equivalent to the dynamic 
yield stress [8]. In general, Y is a functional of the history of deformation 

p 

or any related quantity such as the viscoplastic work W . If the convolution 
form of Stieltjes integral is adopted for such functional 


Y = / e ‘ T/T ° W p ( t - t ) dt 
0 

where t is the time and T 0 stands for the relaxation time constant, and when 
the 3-parameter element model is considered, we get 


Y + aY = Hi W p 


+ H 2 



...(7) 


where a, Hj , and are material constants. 

In order to incorporate the temperature effects into the evolution equations 
(6) and (7), Y is expressed as a function of temperature. Since in the 
proposed constitutive equations Y can be viewed as the equivalent dynamic 
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yield stress, the function may be constructed from experimental data of the 
yield stress versus temperature. A possible form of this function is 


Y 
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Tc - T 
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) 
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where T is a constant, T is a reference temperature, and Y is the value of 

v u u 

the yield stress at T . In addition, from the observation of the variation of 

0 . T 

the plastic flow with temperature, we have £ proportional to exp (-Q/RT ) , 
E ^ exp ( - Q/RT) (9) 


where Q is the activation energy (assumed constant), and R is the universal 
gas constant. From equations (6), (8) and (9) it follows that n must be a 
function of T, 
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The proposed viscoplastic constitutive equations (equations (6) and (7) 
subjected to (8) and (10)) are examined by conducting a series of 
one-dimensional uniaxial numerical calculations. Samples of the results are 
given in Figures 1-4. Figures 1 and 2 display the tensile stress-strain 
results at different strain rates and temperature. Strain rate history effects 
are demonstrated by a jump test in Figure 3. Cyclic test results, Figure 4, 
depict the cyclic hardening effects. Qualitatively speaking, results 
demonstrate the capability of the proposed viscoplastic constitutive equations 
in generating some of the important characteristics of a class of viscoplastic 
materials. A quantitative investigation of the constitutive equations is to be 
conducted experimentally for some viscoplastic materials in a future work. 
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Finite Element Implementation 


Compatibility of the developed coupled thermomechanical equations is 
demonstrated for the case of quasistatic infinitesimal deformation with no 
body force and no heat generation, i.e.. 
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jfo is a scalara constant, p is the density, b stands for the Kronecker 

symbol, a and e are the stress and strain tensors, respectively, 
expressed in vector form, a 1 and £ Vp the corresponding deviatoric stress 
and viscoplastic strain vectors, respectively, e y is the dilatation, and [D] 
is the elastic matrix. 

When adopting the Galerkin finite element method, (11) and (12) become, 
respectively. 


[Kl] 0 + [Cl] t = R + F x ...(13) 

[C2] t + [K2] T = Q + F 2 ...(14) 


where JJ, J, £, and are the nodal displacement, the nodal temperature, the 

nodal force, and the thermal convection load vectors, respectively, F, is a 
vector which accounts for the viscoplastic effects of the balance of momentum 
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equation, and ^ * s a vector which accounts for the mechanical heat 
generation. Also, [Kl], [K2], [Cl], and [C2] are, respectively, the stiffness, 
conductivity, coupling, and consistency matrices. The differential equations 
(13) and (14) can be solved by using the general "9" method in conjunction 
with the fixed point iteration method for the solution of the ensuing 
nonlinear algebraic equations. Results of tensile and compression loading of a 
constrained-ends cylinder for a constant Y can be found in [8]. 

The extension of this work to incorporate hardening and temperature effects 
and solving other practical problems is being undertaken. 
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